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Abstract: 

We review here some general properties of antiferromagnetic Heisenberg spin chains, 
emphasizing and discussing the role of hidden symmetries in the classification of the 
various phases of the models. We present also some recent results that have been 
obtained with a combined use of Conformal Field Theory and of numerical Density 
Matrix Renormalization Group techniques. 

1 Introduction and Summary. 

For quite some time low-dimensional magnetic systems (i.e. (quantum) spins 
on II? and/or 2D lattices) have been considered essentially only as interesting 
models in Statistical Mechanics with no realistic counterpart. It is only in 
recent times that systems that can be considered to a high degree of accuracy 
as assemblies of isolated or almost isolated spin chains and/or of spin ladders (a 
few chains coupled together) have began to be produced and have hence become 
experimentally accessible, thus renewing the interest in their study, which is by 
now one of the most active fields of experimental and theoretical research in 
Condensed Matter Physics. 

In this paper we will discuss only some relevant properties of isolated spin 
chains, referring to the literature [15] for a general review of the properties of 
spin ladders. 

More than one decade ago it was pointed out [20, 31] that integer spin chains 
(more specifically, spin-1 chains, but extensions to different values of the spin 
have also been devised in the literature [43]) possess unexpected and highly non- 
trivial hidden symmetries, whose spontaneous breaking manifests itself through 
the appearance of unusual and highly nonlocal "string" order parameters. The 
string order parameters, together with the more conventional magnetic order 
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parameters, can be used to classify the various phases that the phase diagram 
of one-dimensional magnets can display. 

In the present paper, which is a slightly enlarged version of the talk presented 
by one of us (G.M.) at the XIII — th Conference on "Symmetries in Physics"^ , 
we will concentrate, without pretensions to full generality, on the discussion of 
a few models of antiferromagnetic Heisenberg chains, of their phase diagrams 
and on the role of hidden symmetries in their explanation. The paper is orga- 
nized as follows. In Sect. 2 we review some general facts concerning Heisenberg 
spin chains and discuss how in the continuum limit one can map a "standard" 
(sec below for the terminology) Heisenberg chain onto an effective field theory 
described by a nonlinear sigma-model, and how the presence in the latter of 
a topological term can account for the radically different behaviors of integer 
versus half-odd-integer spin chains. In Sect. 3, concentrating on spin-1 chains, 
we consider the effects of the addition to the " standard" model of biquadratic 
exchange terms and/or of Ising-like as well as of single-ion anisotropics, and how 
the addition of such terms can drive the model away from what is commonly 
called the " Haldane phase" (again, see below for an explanation) towards other 
phases. In this context we will introduce in a more explicit manner the notion 
of hidden symmetries and discuss their role. Sects. 4 and 5 will be devoted to 
the discussion of more recent results that have been obtained by some of us 
[17] with a careful and combined use of analytical (effective actions and Con- 
formal Field Theory) and numerical (Density Matrix Renormalization Group) 
techniques. The final Sect. 6 is devoted to the conclusions and to some general 
comments. 



2 General Features of Spin Chains. 

Let us begin by discussing here what can be considered as the "standard" 

model of an isotropic antiferromagnetic (AFM) Heisenberg chain with nearest- 
neighbor (nn) interactions, which is described by the Hamiltonian: 

N 

n = jY,'Si-'Si+x = n^ + n, (i) 

W_L = {^t^i+i + ^iSt+i) ; '^z = J^'^S^S^j^-^; J_i_ = J^ = J (2) 

i=l 1=1 

where, for each i = 1, ...,N, S'^isa spin operator^: 

, 5f ] = iMijs^^-'S^; a, 7 = x, y, z; 'Sj = t^SiS + 1) (3) 

{S integer or half-odd integer) located at the i — th site of a one-dimensional 
lattice of A'' sites, interacting with its neighbors with an AFM (J > 0) nn 
interaction of strength J. Later on we will consider more general models in 
which J_L ^ Jz will be allowed^. 

^The Conference, organized by Prof. B. Gruber, was held in Schloss Mehrerau in Bregenz 
(Voralberg, Austria), in the period 21-24 July, 2003. 

2 and: = Sf±iSf. 

^J± = 0, in particular, corresponds to the one-dimensional Ising model, a trivially soluble 
classical model. Notice however that an Ising model in a transverse magnetic field becomes 
a genuinely quantum and nontrivial model. 
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It may be useful to define a vector Tij as: "n 



i =: S i/hS, whereby: 





(4) 



Although one is ultimately interested in the thermodynamic (N oo) limit, 
for finite N one can adopt either periodic boundary conditions (PBC's), by 
imposing: 



by which the system is actually considered to "live" on a circle, or open boundary 
conditions {OBC's), where is coupled only to S'2 and S n only to S n-i'^- 
The Hamiltonian of Eq.(l) has an obvious (global) 0(3) symmetry and, for 
P-BC's, it is also invariant under the (discrete) translation group of the lattice. 

In the classical limit (h and S — > 00 with: HS = const.) the spins 
(the r?i's) become (sec Eq.(4)) classical vectors (and lii G the unit 
sphere in K^). The minimum-energy configuration of the spins corresponds 
to: "n j-"rii+i = const. = — 1. Neighboring spins are then aligned antiparallel 
to each other and, in the absence of any external magnetic field, can point in 
a common but otherwise arbitrary direction on the sphere. This is the Neel 
state. Let us remark that, at variance with the /erromagnetic (J < 0) case, 
in which neighboring spins are all aligned parallel, at the quantum level the 
Neel state is not an eigenstate of the Hamiltonian (1). This points to the fact 
that quantum fluctuations will play a much more relevant role in the (quantum) 
antiferromagnetic case than in the ferromagnetic one. 

The classical energy of the Neel state is of course: En = —JN{hS)'^. In this 
state the 0(3) symmetry is spontaneously broken down to 0(2)^, and the state 
exhibits long-range order (LRO). 

Elementary excitations over the Neel state are well-known to be in the form 
of spin waves [38]: coherent deviations of the spins with a dispersion: u){k) oc k 
in the long- wavelength limit {ka <C 1, with a the lattice spacing). Hence, the 
(classical) spectrum of the Hamiltonian (1) is gapless. We would like to stress 
that nothing of what has been said hitherto depends on the value of the spin. 
At the classical level, the spin S ^ can be simply reabsorbed into a redefinition 
of the coupling constant ( J ^ J{hSY) and will contribute only an essentially 
irrelevant and additional multiplicative overall scale factor. 

All this is elementary and well known. Let us turn now to the quantum case^. 
In the early 30's Bcthe [8] and Hulthen [30], employing what has been known 
since as the " Bethe- Ansatz" , were able to show that the quantum S = 1/2 
Heisenberg chain is actually an integrable model. We will not discuss here the 
Bethe- Ansatz in any detail [38], but will only summarize the main features of 
the solution of the S =1/2 model. The (exact) ground state is nondegenerate, 
it exhibits only short-range AFM correlations, but no LRO. Parenthetically, 



*In which case the Hamiltonian sliould be actually rewritten as: Ti. = J^T^ Si - S i^i. 



^Translational symmetry, if present is also broken, as the Neel state is not invariant under 
translations of a lattice spacing as the original Hamiltonian but only of twice the lattice 
spacing. This has important consequences on the location of the Goldstone mode [4, 40, 55] 
in momentum space, that we will not discuss here, however. 

^Or better hS. 

'^Prom now on we will set for simplicity h= 1. 



S i+N = S i^i 



(5) 



JV-l 



i=l 
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this is in agreement with a general, and later, theorem [14]. The (staggered) 
static spin-spin correlation functions: 

g"{i-j) = {-iy'-^^{SrSf); a = x,y,z (6) 

where (...) stands for the expectation value in the ground state, are all equal 
and decay algebraically to zero at large distances. We recall here that genuine 
LRO would imply (we omit here the index a): 

lim g{i-j) = ON^0 (7) 

this defining the Neel order parameter On (actually the square of the equi- 
librium staggered (i.e sublattice) magnetization). On the other extreme, an 
exponential decay of the correlations of the form, say: 

g{i-j) ~ exp{-|i-j|/^}P(H-j|) (8) 

with P{.) some inverse power of \i — j\ would imply a finite correlation length 
^ and a mass gap (or, better, a spin gap) A in the excitation spectrum roughly 
given by: A oc c^~^,with c a typical spin- wave velocity. Algebraic decay of 
correlations (formally corresponding to ^ — > oo) implies then that the system 
is gapless. Summarizing, the main features of the S ^ 1/2 Heiscnberg AFM 
chain are that it has a (quantum) disordered ground state, with only short- 
range AFM correlations, and that it is gapless. It is therefore a (actually the 
first) prototype of a (quantum disordered and) quantum, critical system [48]. It 
can be said then that, as compared with the classical limit, the system remains 
gapless but quantum fluctuations destroy LRO. 

About thirty years later Lieb, Schutz and Mattis [36] (LSM) proved an 
important theorem stating that an S = 1/2 chain has either a degenerate 
ground state or is gapless. No surprise that the Bethe solution obeys the Lieb- 
Sehutz-Mattis theorem, which is however of much wider reach, as it can cover 
models that are more general than the "standard" nn chain, such as, e.g., the 
Majumdar- Ghosh [37] model, another integrable model that we will nor discuss 
here, though. The results of LSM were extended later on by other authors 
[3] beyond S = 1/2 to cover all the half-odd- integer values of the spin. One 
can then take as rigorously proven that (at T = 0) isotropic half-odd-integer 
Heiscnberg chains (with constant nn interactions) are all quantum, disordered 
and quantum critical (i.e. gapless). This result was thought for quite some 
time to be "generic", i.e. valid for chains of arbitrary spin until, in the early 
80's, Haldane [28] put forward what has become known since as Ha,ldane's 
conjecture'^ , according to which half-odd-integer spin chains should be quantum 
disordered and gapless but integer spin chains should instead exhibit a spin gap 
and an exponential decay of correlations. This implied that, contrary to what 
happens in the classical limit, the physical behavior of spin chains should be a 
highly discontinuous function of the value of the spin. 

Completely rigorous proofs of (the second part of) Haldane's conjecture are 
still lacking. However, strong support to it comes from the analysis of the 
continuum limit of the Heiscnberg chain, which we will briefly describe now, 
referring to the existing literature [1, 6, 23] for more details. 
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The canonical partition function for the Hamiltonian of Eq.(l) at tempera- 
ture T = {kBl3)~^ (with ks the Boltzmann constant): 



3 = rr{exp[-/3H]} 



(9) 



can be written as a spin coherent-state path-integral [32], whereby the spin 
variables get replaced, inside the path-integral, by classical variables according 
to: 

5i ^ SQi (10) 



with fli a classical unit vector: 



1. The next (and perhaps the most 



important) step in Haldane's analysis is the parametrization of the Oj's as^: 



h. 

s 



+ 



(11) 



with: \ni\ = 1 and: rii- I i = 0. The n^'s are assumed to be slowly- varying 
(on the scale of the lattice spacing). In this way, capitalizing, so-to-speak, on 
the information gained from the Bethe-Ansatz solution of the 5 = 1/2 model, 
they incorporate the information that the system still retains some short-range 
AFM ordering, which would be global only for rii = const, (and ^ j = 0). The 
I i's can be shown [1] to be the (local) generators of angular momentum. In 
the semiclassical (large S) limit, an expansion of the action in the path-integral 
up to lowest (second) order in the I i's is justified. Taking then the continuum 
limit together with a gradient expansion, and integrating out the I j's, one ends 
up with the following expression for the partition function: 



3 = y" [Vn] S (n^ _ 1) exp{-5B - iSe} 



(12) 



where [Dn] stands for the functional measure and the 5 inside the integral is a 
functional 6. The first term in the action is given by: 



Se 



dx I dr 







29 



-\drn\'^ + c\dxn\'^ 



; n = n{x, r) 



(13) 



where L{= A'' x lattice spacing) is the length of the chain, g = 2/S is the coupling 

constant and: c = 2JS is the spin-wave velocity. This is simply the Euclidean 
action of an 0(3) nonhnear sigma model [6, 23, 59] {NLaM). The second term 
is the integral of a Berry phase [50] , and is given by: 



L (i 



Sb = 



4n 



dx J drn- {drnxOxfi) 



(14) 







with: 6 = 2ttS. The coefficient of 6 is easily recognized to be the Pontrjagin 
index [11, 41, 44], or winding number, of the map: 



b2 

comp 



(15) 



'This is what is known in the Uterature as "Haldane's mapping" 



5 



from spacctime compactificd to a sphere and the two-sphere where n takes 
values, and it is an integer: Sb is therefore a topological term, and: Sb = 2TmS, 
n e Z. Therefore, exp{— i^s} = 1 for integer S {9 = mod(27r)), but: 
exp{— zS'b} = (—1)" {6 = TT mod (27r)) if the spin is half-odd integer. This will 
generate interference between the different topological sectors, and it is the at 
the heart of the different behaviors of the two types of chains. 

The pure {6 ~ in our case) (1 + 1) 0(3) NLaM is a completely inte- 
grable model [60]. It has a unique ground state, and the excitation spectrum 
is exhausted by a degenerate triplet of massive excitations that are separated 
from the ground state by a finite gap. On the contrary, the 6 = tt model was 
shown [45] to be gapless. Therefore, Haldane's conjecture is fully confirmed by 
the analysis of the continuum limit of the Heisenberg model. 

We would like only to mention in passing that quite a similar behavior occurs 
in spin ladders [15, 19, 46], namely even-legged ladders are gapped, while odd- 
legged ladders are gapped for integer spin and gapless for half-odd-integer spin. 
This ^' even-odd^' effect has been shown [19, 51] to have the same topological 
origin an in single chains. 

How do these results compare with the gaplessness (irrespective of the value 
of S) of the 5*^00 classical limit? The answer resides in the dependence of 
the spin gap on S. Already at the mean-field level, but more accurately from 
large- A'' expansions and/or renormalization group analyses [42], it turns out 
that the spin gap A behaves as: A cx exp{— ttS'} for large S ^. Hence, integer- 
spin chains become exponentially gapless for large S, and the classical limit is 
recovered correctly. 

3 More general Models. Hidden Symmetries 
and String Order Parameters. 

In view of what has been said up to now, the second part of Haldane's conjecture 
is by far the most intriguing part of it. Therefore integer spin chains are the 
most interesting ones, and we will concentrate from now on on iS = 1 chains. 

What has been called in the previous Section the " standard" AFM Heisen- 
berg model is actually a member of at least two larger families of models that 
we will illustrate briefly here. The first class of models, that we will call 
"6— models", includes a biquadratic term in the spins, and is described, setting 
J = 1, by the Hamiltonian: 

H = cose pi- + sine pi-^.+i) [ (16) 

i=l ' 

with = corresponding of course to the " standard" model. Most of the phase 
diagram has been obtained [31] numerically, except for the points at 9 = ±7r/4, 
that correspond to integrable models. The point e = 7r/4 is the Sutherland 
model [53], while 9 = —n/A is the integrable model [9, 54] of Babujian and 
Takhatajan^''. Both models are gapless, while the entire region — 7r/4 < 9 < 7r/4 
is known (numerically, again) to be gapfull. This whole region has been called 

oc exp{— irnfS} for spin ladders [19], where n; is the number of legs of the ladder. 
l^This point is also known familiarly as the "Armenian point". 
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the " Haldane phase^^ . It includes a particularly interesting point that has been 
studied extensively by Affleck, Kennedy, Lieb and Tasaki [2] {AKLT), namely 
6 = 6*, with: tan^* = 1/3. The corresponding Hamiltonian (omitting an 
irrelevant overall numerical factor) is given by: 

'Haklt = I S j- S i+i + (s j- S i+i] > (17) 

This model is not completely integrable, but the ground state is known, 
it is unique in the thermodynamic limit and can be exhibited explicitly. The 
ultimate reason for this is that, apart from numerical constants, the i — th term 
in curly brackets is just: 



S i- S i+i + - (^S i- S 



2 



(18) 



where P2{i, z + 1) is the projector [39] onto the state of total spin Stot — 2 of the 
pair of 5 = 1 spins located at sites i and i + Therefore, the ground state of 
Ti-AK LT must lie in the sector of the Hilbert space that is annihilated by all the 
projectors. It was shown by AKLT that the exact ground state (also called the 
"Valence-Bond-Solid" {VBS) state) can be constructed as a linear superposition 
of states "I>CT that have the following characteristics. Let: a = {(Ti, ...cttv} be a 
given spin configuration, with: Ui = 0, ±1, t = 1, A^. Then, is such that: 
i) = CTj^CT and moreover: ii) If a given spin is, say, +1, then the next 

nozero spin must be —1, and viceversa. Typical such states correspond therefore 
to spin configurations of the form: 

CT = {-hOO-0-F-000+...} (19) 

In other words, "up" and "down" spins do alternate in <I>cr, but their spatial 
distribution is completely random, as an arbitrary number of zeroes can be 
inserted between any two nonzero spins. So, if a given spin in nonzero, we can 
predict what the value of the next nonzero spin will be, but not where it will 
be located. There is therefore no long-range (Neel) order in any conventional 
sense in the VBS ground state, but a sort of "Liquid Neel Order" [LNO). 
Conventional Neel order would be characterized by a nonvanishing of (at least 
one of) the Neel order parameters: 

0%= lim (-l)l*-^l ; a = x,y,z (20) 

\l-j\^<X 

In the VBS state and (numerically) in the whole of the Haldane phase one finds 
instead [2, 31]: = 0,a = x,y, z, and this is consistent with the absence of a 
"rigid" Neel order. 

There remains however what we have called the "liquid" Neel order, and it 
has been argued convincingly in the literature [20, 31] that this is connected 
with the nonvanishing of a novel class of order parameters that we will discuss 
now briefly. Let us begin by defining the string correlation functions as: 



g^{n) =:-(5o«exp 



■ n-1 



S^); a = x,y,z; n>0 (21) 
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These are similar to the standard two-point correlation functions: 



e«(n)=:(-l)"(So«5«) (22) 

whose asymptotic (n — > oo) limit yields the Neel order parameter (s), except 
that a string of exponentials of intermediate spins has been inserted between 
the leftmost and the rightmost spins. 

The string order parameters (SOP's) are then defined as: 

Og= limgg(n) (23) 

It turns out [2, 25] that the string correlation functions are strictly constant 
in the AKLT ground state, namely: 

ag(n) = canst. = = ^ (24) 

The ground-state spin-spin correlation functions have also been evaluated 
exactly for the VBS state [2] , and they turn out to be given by: 

^"W = ^f^V (25) 



3 V3, 

In other words: Q"{n) oc exp {— ti/^aklt}, where the correlation length ^aklt 
is given, in units of the lattice spacing, by: 

Uklt = 0.91 (26) 

log 3 

less than unity in units of the lattice spacing, implying a rather large spin gap. 

So far for the ground state of the AKLT model. String and ordinary cor- 
relation functions as well as Neel and string order parameters have also been 
evaluated (numerically away from the AKLT point) for other points of the Hal- 
dane phase [25]. For example, at the Hciscnbcrg point, exact diagonalization 
methods^ ^ have shown that the string correlation functions are not strictly con- 
stant, but still decay exponentially to a value of the string order parameter that 
is somewhat smaller (Og ~ 0.36...) than the AKLT value (Og = 4/9 ~ 0.44...) 
but still nonzero. The spin correlation length was also found [25] to be slightly 
larger than the AKLT value, but still finite. So, there is convincing evidence 
that the entire Haldanc phase is characterized by vanishing Neel order param- 
eters but by nonzero SOP's. There is also convincing numerical evidence [25] 
that the string order parameters vanish at the integrable boundaries of the 
Haldanc phase (i.e. for 9 = ±7r/4). 

That the nonvanishing of the SOP's is connected to the breaking of a sym- 
metry, and hence to the onset of an ordering that is not apparent in the original 
Hamiltonian was clarified in a seminal paper by Kennedy and Tasaki [31] {KT). 
With reference to a given configuration {u}, and defining N{a) as the number 
of odd sites at which the spins are zero, one defines a new configuration {a} via: 



Ui = exp 



i-l 



(^i (27) 



^'^With the Lanczos method and for chains with up to no more than 14 sites. 
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and then a unitary operator U via: 



= (-l)^('^) $^ (28) 

In a nutshell, the action of U amounts to leaving the first nonzero spin unchanged 
and to flipping every other nonzero spin proceeding to the right of the chain. 
For example: 

{+ + - +00 + - + +} {+ - - -00 + + + -} , . 
{0 + -00 + 00 - +00-} {0 + +00 + 00 + +00+} ^ ^ 

and so on. It is obvious that J7 is a unitary-'^^. What is less obvious is that 
the unitary transformation is a nonlocal one, in the sense that U cannot be 
written as a product of unitary operators acting at each single site. This has 
the important consequence that symmetries that are local (in the above sense) 
for the Hamiltonian 7i will of course remain symmetries of the transformed 
Hamiltonian H (as U is unitary) but need not survive as local symmetries of 
H. Specifically, the symmetry group of H is SU{2), that includes a discrete 
Z2 X Z2 subgroup of rotations of n around the coordinate axes. Explicitly, the 
transformed Hamiltonian has the form [31]: 

'H = J2 {cos Ohi + sin (hi) } (30) 

i 

where: 

hi = -5f 5f+i + 5?' exp {in {Sf + SU,) } Sf^, - SfSf^, (31) 

and it evident that Z2 x Z2 is the only local surviving symmetry group of the 
transformed Hamiltonian Ti. Even more important is how the string order 
parameters transform. The result is [31]: 

0% in) = o%„, (n) (32) 

where: ^ 

0%rro{n) = {StSf)\^ (33) 

and the r.h.s stands here for an average taken w.r.t. the ground state of the 
transformed Hamiltonian. The transformed order parameter is now a ferromag- 
netic order parameter. Therefore: Og {H) ^ =J> Cje^^o (^) 7^ 0' ^'^'^ ^^^^ 
implies the onset of a spontaneous ferromagnetic polarization in the a — th direc- 
tion in the ground state of Ti. This in turns entails a partial (if O"^^^^ (l-tj ^ 
for just one value of a) or total (if this happens in more than one direction) 
spontaneous breaking of the discrete Z2 x Z2 symmetry. It is known [4, 55] that 
spontaneous breaking of a continuous symmetry is accompanied by massless ex- 
citations (the Goldstone modes) , while breaking of a discrete symmetry usually 
implies the opening of a gap (the most conspicuous and familiar example being 
the 2D Ising model). Therefore, KT were led to consider the spontaneous 
breaking of the Z2 x Z2 symmetry as the origin of the Haldane gap. 

^■^ Notice also that: N (ct) = AT (ct), as zero spins are mapped into zero spins. 
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One has however to be a bit careful on this point. It appears to be true 
that spontaneous (partial or total) breaking of the Z2XZ2 symmetry implies 
the generation of a spin gap. But: 

i) The converse need not be true. We will see that there are spin models 
that exhibit gapped phases^^ while retaining the full Z2 x Z2 symmetry, and: 

a) The mere nonvanishing of (one or more) string order parameters is not 
enough to fully determine in which (gapped) phase the system is. It is the full set 
of order parameters, both string and Neel, that allows for a full characterization 
of the various phases. In particular, the Haldane phase is fully characterized 
by the vanishing of all the Neel order parameters and by all the three string 
parameters being nonzero. 

We turn now to a different class of models, the so-called "A — D" family of 
models^^. They arc described by the family of Hamiltonians (parametrized by 
two real parameters, A and D): 

N 

'H^{s^S^+,+S!Sf^,+\S^S!+,+D{Slf] (34) 

The "standard" (isotropic) AFM Hciscnborg model corresponds of course to 
A = 1 and D = 0. A = — 1 (and D = 0) can be easily shown^^ to correspond to 
a (isotropic) ferromagnetic Heisenberg model. |A| 7^ 1 introduces an "Ising-like" 
anisotropy, while D introduces what is called "single- ion" anisotropy. 

The model can be solved exactly for S = 1/2 [33], but no exact solutions 
are available for integer spin. There are obvious asymptotic limits when either 
A (resp. D) is large and D (resp. A) not too large, so that the "A-term" (resp. 
"f-term") can be considered as a zeroth-order Hamiltonian and the rest as a 
perturbation: 

i) |A| 1. The reference ground state is either a Neel AFM state (A > 0) 
or a ferromagnetic (A < 0) state. 

a) \D\ ^ 1. For D > (the so-called "large- 1?" phase) the reference state 
becomes a planar state with Sf = for all i's, while for _D < the reference 
state is a state where S^ = is excluded, hence a state where the S = 1 spins 
become effectively two-level systems, and a detailed map of the model into an 
effective spin-1/2 model [17, 47] can be successfully performed. For |A| ^ 1 
and D the symmetry group of the Hamiltonian is 0(2)xZ2 (the Z2 factor 
corresponding to a reflection in the x — y plane: 5? —Sf). 

Apart from these limiting cases, the model has been studied analytically [49] 
as well as numerically [10, 13, 26, 52] , and the corresponding phase diagram is 
displayed in Fig.l. 

The various sectors of the phase diagram can be characterized as follows 
[13, 17, 22]: 

i) the Haldane phase. The ground state is unique with total component 
S^g^ = of the spin. The order parameters are: = 0, Og 0,a = x,y, z. 
The (Haldane) gaps are in different spin channels according to the sign of £>, 
but nonzero in any case. The isotropic, Heisenberg point A = 1, £> = is in this 

^'^The so-called "large-D" phases of the "A — £>" model to be discussed immediately below. 
^*This is the class of models on which the Bologna group is currently working. 
'^^By performing a rotation of tt around the z-axSs on one of the two sublattices (i.e. on 
every other site). 
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Figure 1: Phase diagram of the X — D spin-1 Hamiltonian of Eq. (34). The 
indicated regions are explained in the text. 



phase, and Ues on a Une separating the two subphases, that are denoted as HI 
and H2 in the literature [10]. 

a) The Neel phase. The ground state is doubly degenerate, and the order 
parameters are: = Og = for a = x, y, but: O^, Og ^ 0. 

Hi) The large-D phase.The ground state is unique, it is gapped, but here:Og = 
O% = Va. 

iv) The two XY phases. These are both gapless phases. They are dis- 
tinguished by the nature of the low-lying spin excitations (spin-1 in the XYl 
phase, spin-2 in the XY2 phase). 

v) The ferromagnetic phase. The ground state is doubly degenerate, with 

maximal magnetization: ^f^^ = ±A'^, and the phase is gapped. In this case 
it is the /erromagnetic order parameter that is nonvanishing, and actually [16]: 



^ ferro — ^ (^^^ Other two being zero). Also: Og{j,k) = (— )•' ^, while: 
C^'^ = 0. 



Anticipating some of the numerical results of Sect. 5, we give below, in Figs. 2 
and 3, some examples [IG] of the behavior of the various correlators and order 
paranictc;rs as functions of the parameters of the model. 

Concerning the nature of the transitions between the various phases [10, 13, 
17], both the Haldane-large-i? and the Haldane-Neel transition lines are critical 
(gapless) lines. The two critical lines merge at a tricritical point (at I? ~ A ~ 3), 
above which the Haldane phase disappears and the transition (a large- iP-Neel 
transition,now) is first-order. The XF-ferromagnetic transition is instead a 
first-order one, as well as the large-D-ferromagnetic transition. Finally, The 
Haldane-Xy transition is considered [13] to be of the Berezinskii-Kosterlitz- 
Thouless [7, 34] [BKT) type, as well as the XY-\axge-D transition. 
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Figure 2: Ordinary and string correlation functions in the Haldane phase: (a) 
g^{k), (b) 5(-)'=(S'c^5'^:), (c) a|(A;) and (d) g%{k). Selected values of the param- 
eters are (i? = 0.5,A=l). Note that with this choice the transverse correlation 
length is appreciably larger than the longitudinal one. The data have been ob- 
tained with finite-size DMRG on a chain of L = 100 spins {S = 1) with PBC 
and M = 216 states (Sect. 5 for details). 



The "A — £)" model has also been studied by KT. Applying the same 
nonlocal unitary transformation that was discussed previously, they showed that 
the transformed Hamiltonian, whose explicit form we will not give here, is still 
given in terms of the operators hi (see Eq.(31)), and retains therefore Z2XZ2 
as the only local symmetry,just as in the case of the Hamiltonian of Eq. (16). 
Therefore, the same conclusions as before apply concerning the connection of 
the nonvanishing of the string order parameters with the spontaneous breaking 
of the Z2 X Z2 symmetry. 

In the present paper we will address mainly to the detailed nature of the 
Haldane-largc-Z) and Haldane-Neel critical transition lines. It is known that 
the (large distance) critical behavior of one-dimensional quantum systems is 
well described by Conformal Field Theory [12, 21, 24, 27, 33] {CFT). In the 
next Section we will report on a proposal of an effective CFT for the "A — D" 
model on the Haldane-large-D critical line. This allows for the prediction of the 
operator content of the theory, and hence also for the prediction of the structure 
of the conformal tower of excited states above the ground state. To confirm the 
predictions, we will report also on extended numerical analyses, whose details 
will be reported elsewhere [18], that fully confirm the theoretical predictions. 

4 Conformal Field Theory and Effective Actions. 

Let us b(;gin by recalling some basic results and examples of CFT that will be 
used in the forthcoming analysis of the critical properties of the spin-1 X — D 
chain. 

It is well known [21, 27] that critical properties of two-dimensional systems 
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Figure 3: Order parameters relevant to the Nocl-Haldanc-largc D transitions 
plotted versus the anisotropy coefficient D of Eq. (34) fixing A = 1: (a) 
defined in Eq. (20); (b) 0% and (c) 0% defined in Eq. (23). The asymptotic 
values are extrapolated using an algebraic best-fit fimction Oqo + C /\i — on 
the DMRG data (same choices as in Fig.2). Near D ~ —0.3 Oj^ and Og do not 
vanish in the same point due to finite-size effects and to a moderate number of 
DMRG states. 



are completely classified by CFT's: since in 2D the conformal group is infinite 
dimensional, the Hilbcrt space of a conformally invariant theory can be com- 
pletely understood in terms of the irreducible representations of its algebra, the 
Virasoro algebra. We recall that the latter has an infinite number of genera- 
tors, denoted with L„, Z„ {n E Z) for its holomorphic and antiholomorphic part 
respectively, satisfying the commutation relations: 

[Ln, Lm] = (m - n)Lm+n + ^("^^ " rn)Sm+n,o (35) 

and similarly for the Z„. The constant c is called the central charge of the alge- 
bra or the conformal anomaly. Since we are interesting in a comparison between 
theoretical predictions and numerical data, which are performed on a finite lat- 
tice, we will consider a CFT defined on a cylinder with spatial dimension of 
finite length L. In this case [12, 21], the energy and the momentum operator 
are represented respectively by: 

H = ^(^o + Lo-^) (36) 

P = ^(Lo-Lo) (37) 

In order for H to be bounded from below, we must restrict our attention to 
highest weight representations of the Virasoro algebra, for which there exists a 
highest weight (or primary) state |A, A) satisfying: 

Lo|A, A) = A|A, A) , L„|A, A) = for n > (38) 
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and analogous relations with rcpcct to the X„ generators. 

Each of these representations is thus identified by the values of the central 
charge c and of the couple (A, A) (the conformal dimensions). They fix both 
the energy and the momentum of the primary state | A, A), according to: 

Ka = + (39) 

^A,A = (40) 

Notice that, in a finite geometry (with PBC), the vacuum state, corresponding 
to A = A = 0, has a non zero energy (Casimir effect): 

^0 = - — (41) 

Also, the two-point correlation function of the operator creating a given primary 
state out of the vacuum (| A, A) — a|0)) has an algebraic decay whose critical 
exponents are determined by the values of the conformal dimensions (A, A): one 
has [21, 27] 

{0^,a{z, z)O^M^, 0)) a (42) 

Finally, from the primary state |A,A) one can obtain all excited (or sec- 
ondary) states by applying strings of powers of L„ with n < 0. It is easy to 
see that, if m,n < 0, the commutation relations (35) imply : Lo(-t'm)'' |A, A) = 
(A + mj)|A,A), Lo{Ln)''\A,A) = (A + nfc)|A,A), so that the secondary states 
have energies and momenta: 

E^S-Elc = ^(A + A + r + f), (43) 

<A = ^(A-A + r-f) (44) 

with r, f G N and a degeneracy that can be explicitely calculated for each 
representation. It may happen that some of these states have null norms. In this 
case the true (non-degenerate) Hilbert space of states is obtained after projecting 
out these null vectors, which therefore do not contribute to the operator content 
of the corresponding CFT. The quantity in brackets in the right hand side of 
Eq. (43) yields the coefficients with which the energy of the corresponding 
state scales to zero in the thermodynamic limit. It is therefore called "scaling 
dimension" and will be denoted by (i^'^ in the sequel. 

Let us examine some examples. We will consider only unitary theories, 
corresponding [21] to the following two sets of values of the central charge c: 

c = 1-^,^ = 3,4,...; (45) 

c > 1 . (46) 

The first set of values corresponds to the so called minimal models [21, 27], 
whose primary states are of finite number. Their conformal dimensions are 
given by the formula: 

A,„A,,= [(P+l>-H'-l ,i<s<r<p-l,r,s€Z (47) 
4p{p + 1) 
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Theories with c > 1 have instead an infinite number of primary states. 

The simplest case of a CFT corresponds to c = 1/2 (p = 3 in Eq. (45)) and 
describes the univerahty class of the two-dimensional Ising model. According to 
(47), there are only three primary operators: the identity I corresponding to the 
vacuum, (A, A)i = (0, 0), the Ising spin a with (A, A)o- = (1/16, 1/16) and the 
energy density e with (A, A)^ = (1/2, 1/2). Notice that the spin-spin correlator 
(CT(x)cr(O)) decays with a critical exponent 77^ = AA^ = 0.25. In Table 1 we list 
the lowest conformal (primary and secondary) states, together with their scaling 
dimensions and momenta. As explained in the next section, a comparison with 
the numerical data given in the last column will allow us to conclude that the 
Haldane-Neel critical transition line is indeed of the Ising type. 



Table 1: Columns 1-4 show the conformal dimensions (A, A),(r, r), the scaling di- 
mensions d^'^' and the momenta P^'^ of the lowest conformal states in the c = 1/2 
minimal model. The numerical results in the last column are explained in Sect. 5. 
Notice that the states with A = A = 0, (r, f) = (1, 0), (0, 1) do not appear since they 
correspond to null vectors. 



(A, A) , (r,f) 


Ar,r) 

"a,a 


p(r,r) 




(0,0); (0,0) 










(1/16,1/16): (0,0) 


1/8 





0.1250 ±0.0004 


(1/2, 1/2); (0,0) 


1 





0.962 ±0.001 


(1/16,1/16); (1,0), (0,1) 


9/8 


±27r/L 


1.0959 ±0.0008 








1.100 ±0.003 


(1/2,1/2); (1,0), (0,1) 


2 


±27r/L 


1.87 ±0.02 








1.87 ±0.02 


(0,0); (2,0), (0,2) 


2 


±47r/L 


1.904 ±0.004 








1.86 ±0.01 



We discuss now briefly the c = 1 case, which exhibits a much richer structure. 
It corresponds to the field theory of a free compactified bosonic field, i.e. to a 
Gaussian model with Lagrangian: 



(48) 



where 9 represents an angular variable spanning a circle of a given radius R and 
the constant v, which has the dimension of a velocity, is called spin velocity. If we 
assume for Q, and hence for its dual field periodic boundary conditions, the 
Hilbert space of the theory splits into a direct sum of distinct topological sectors 
labeled by the winding numbers n, m G Z of the fields 9 and $ respectively. 
The primary fields are then vertex operators of the form [21, 27] 

Vmn = exp (^WATrKn^ + i^-K/KmO^ (49) 

whose scaling dimensions are given by 

dmn={^+n'K\, K=^ (50) 



^^If we decompose the field © in its iiolomorpiiic and antiholomorphic part, Q{z,z) 
&h(z) + Qah(z), the dual field is defined as $ = — Qahi^)- 
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Notice that the latter depend exphcitely on the radius of compactification. Thus 
we obtain a different c = 1 theory for each value of R, i.e. of K. For example, 
K = 1 corresponds (via fermionization [21, 27]) to a ID model of free Dirac 
{FD) fcrmions. The K = 1/2 point is said to be self-dual (SD) since it is 
invariant under the duality transformation 6 <^ m <^ n, while the point 
K = 2 corresponds to the BKT critical theory. 

We remark also that the energy operator (90)^ has conformal dimension 2 
for any value of R and hence it is always marginal. The effect of adding it to 
the Lagrangian (48) results only in a change of the coupling constant in front, 
which, in turn, can be absorbed into a rescaling of the radius of compactification 
of 9. Thus we generate a continuous line of inequivalent critical c = 1 theories, 
corresponding to different values of K. 

It is well known [27] that the Gaussian model (48) describes the contin- 
uum limit of the spin 1/2 XXZ chain with anisotropy parameter A, as long as 
— 1 < A < 1. Prom the exact Bethe-ansatz results, one can show [27] that the 
interesting cases A = —1, 0, 1 corrrcspond to the SD^ FD and BKT points of 
the bosonic theory, respectively. We would like to show now, that the Gaussian 
model (48) describes also the critical properties of the spin-1 \ — D Hamiltonian 
(34) on the Haldane-largc-D transition line. In doing so, we will also establish a 
relationship between the coupling constants D, A of the discrete model and those 
of the continuum theory, namely the spin-wave velocity v and the compactifica- 
tion radius. This will allow us to make quantitative theoretical predictions to 
be compared, in next section, to the numerical results. 

In the spirit of Haldane'a mapping, we start from a classical solution, which 
for D > A — 1, is a planar state where the unit vectors ^jir) that represent our 
spins {S j Silj{T) , see Sect. 2) are Neel ordered in the xy-pl&ne: ilj{T) = 
(cos((^o + JTr), sin(6'o + jtt), 0). Hence we make the Haldane-like ansatz: 



hj{r) = (-l)%,(r)y (51) 

where hj{T) = e'^'*^^' e 0(2)^^^, z is the unitary vector (0,0,1), and the fluc- 
tuation field Ij is supposed to be small. Thus, as for the isotropic case, it is 
possible to obtain an effective Lagrangian that describes the low-energy physics 
of the Hamiltonian (34) in the continuum limit. Carrying out this calculation 
as explained in Sect. 2, one obtains in this case a Gaussian model (48), where 
now O = 0/^ and 

g=^y/2{l+D + X)- v = s^/2{l + D + X) (52) 

In other words, we have a free theory for a bosonic field O, which is compactified 
along a circle of radius Thus, the operator content of the theory can be 

read from Eq. (49): the list of primary operators is exhausted by the vertex 
oprators Vmn whose scaling dimensions are given by Eq. (50), with K = ir/g. 

In addition, the scaling dimensions (50) fix also the (non universal) critical 
exponents of the correlation functions. For instance it is easy to see that the 
transverse spin-spin correlator should decay according to: 

{S+{0)S-{x)) « (e*«(o)g-ie(x)jj ^ |^|-r, ^^^^ ^ ^ ^dw = q/^tv . (53) 
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5 The Density Matrix Renormalization Group 
and Spin Chains. 

The code that we have used for density matrix renormalization group {DMRG) 

calculations follows rather closely the algorithms reported in White's seminal 
papers [56, 57], with the following points to be mentioned: 

• The superblock geometry was chosen to be [B* • |B^gf»] with PBC, where 
B^^f is the (left ^ right) reflected of block B* with ,s' sites. The rationale for 
adopting this configuration is that, being effectively on a ring, the two blocks 
are always separated by a single site, for which the operators are small matrices 
that are treated exactly (no truncation) [57]. In this way we expect a better 
precision in the correlation functions calculated fixing one of the two point on 
these sites and moving the other one along the block. Moreover, whenever the 
system has an underlying antiferromagnetic structure (typically when a stag- 
gered field is switched on), this geometry seems to be the one that preserve it 
at best, both for even and odd values of s. 

• We used the finite- system, algorithm, with three iterations. This prescription 
should ensure the virtual elimination of the so-called environment error [35], 
which is expected to dominate in the very first iterations for L < L*{m) (see 
below). Normally the correlations are computed at the end of the third itera- 
tion, once that the best approximation of the ground state is available. This has 
the advantage of using less memory during the finite-size iterations but requires 
the storage of all the matrices needed to represent, on the reduced basis of the 
last step, the operators entering the correlation functions of interest. At the 
moment, disk storage is the ultimate factor that limits the size of the systems 
that we are able to treat. 

• We always exploit the conservation of S^^^. With the exception of the fer- 
romagnetic phase, that we do not address now, the ground state(s) is (are) at 
^tot — [10]. In order to maximize their accuracy, the correlations are cal- 
culated targeting only the lowest-energy state within this sector. However, in 
order to analyze the energy spectrum, we had to target also the lowest-energy 
states in the other sectors IS'foJ = 1,2, .. . and/or a few excited states within 
the S^g^ = sector, depending on the phase of interest. On the one hand, this 
requires a modification of the basic Lanczos method to go beyond the lowest 
eigenvalue of the superblock Hamiltonian. On the other hand, once the A^t 
eigenvalues of interest are found, one can build the block density matrix as the 
average (mixture) of the matrices associated with the corresponding eigen- 
vectors. At present we are not aware of any specific "recipe" other than that of 
equal weights. 

Going back to the modified Lanczos routine, our DMRG code implements the 

so-called Thick Restart algorithm of Wu and Simon [58]. Once S^^^ is fixed, in a 
given run we want to determine simultaneously the first A't levels \S^g^]h) with 
b=0,l,2,. . . , A/t - 1 (the ground state being identified by {S^^t = 0,b= 0)). Then, 
as in the conventional Lanczos scheme, we have to push the iteration until the 
norms of the residual vectors and/or the differences of the energies in consecutive 
steps are smaller than prescribed tolerances (10~^ — 10~^^ in our calculations). 
The delicate point to keep under control is that, once the lowest state \S^g^,Q) 
is found, if we keep iterating searching for higher levels the orthogonality of the 
basis may be lost, just because the eigenvectors corresponding to these levels 
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tend to overlap again with the vector |5foj;0). As a result, the procedure is 
coniputationally more demanding to the extent that one has to re-orthogonalize 
the basis from time to time. Typically, we have seen that this part takes a 
10-20% of the total time spent in each call to the Lanczos routine. We have also 
observed that if this re-orthogonalization is not performed, one of the undesired 
effects is that the excited doublets (generally due to momentum degeneracy) are 
not correctly computed. More specifically, it seems that while the two energy 
values are nearly the same in the asymmetric stages of the iterations, when 
the superblock geometry becomes symmetric (s = s' in the notations of the 
preceding point) the double degeneracy is suddenly lost and only one of the two 
states appears in the numerical spectrum. 

So far for the specific algorithm. Now, the crucial point to consider in accu- 
rate DMRG calculations is the choice of M, that is, the number of optimized 
states. White argued [57] that the convergence of the ground state energy is 
almost exponential in M with a step-like behaviour, probably related to the 
successive inclusion of more and more complete spin sectors. Unfortunately, 
the effective accuracy gets poorer when we deal with energy differences and 
correlation functions, for which little is known about convergence. It must be 
told, however, that despite its name the DMRG performs somehow better for 
systems with a definite gap rather than for gapless (critical) ones. We refer to 
the papers by Andersson, Boman and Ostlund [5] and by Legeza and Fath [35] 
where, for different systems and in terms of different observables, the following 
common feature emerges: Even if the quantum system is rigorously critical in 
the limit i — > oo, the DMRG truncation introduces a spurious length, L*{m), 
which, as expected, diverges as M is increased. (Our analysis of the accuracy 
of the energy levels in some selected points of the X — D chain near criticality 
leads to a similar conclusion [18]). Hence, even if we are technically able to 
deal with sizes L > L*{m) (at a given M), as far as criticality is concerned we 
cannot rely completely on the DMRG data because the system experiences an 
effective length which should be absent in the critical regime. 

Therefore, our strategy can be summarised as follows: We fix a rather high 
value of M such that the trustable values of L are sufficiently large to see 
the scaling limit of CFT, but not too large as compared to L*{m). In other 
words, even in the study of (supposed) critical systems, we prefer to exploit 
the computing resources to include as many DMRG states as possible, and 
to refine the calculations with finite-size iterations, rather than trying to take 
naively the limit L oo. In addition, to judge whether M is sufficiently large 
or not we checked the properties of traslational and reflectional invariance that 
the correlation functions should have ^"^ . To be specific, if G{0,k) is a certain 
correlation function computed starting at j = 0, we have always increased M 
(at the expenses of L) until the bound \Q{£,£±k) — G{0, k)\/G{0, k) < 0.05 was 
met for k varying from to £ = L/2, possibly with the exception of the ranges 
where G{0, k) is below numerical uncertainties (10~®, say). 

The quality of the numerical analysis of the critical properties depends heav- 
ily on the location of the critical points of interest. As far as the transitions 
from the Haldane phase are concerned, it is convenient to fix some representa- 
tive values of A and let D vary across the phase boundaries. This preliminary 

'^'^In [18] it is shown that, in general, Gg'^{j,k) behaves nontrivially under j <-> k, due to 
the fact that the ground state is not necessarily in the S^^f = sector. 
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task of finding Dc{X) turns out to be crucial for subsequent calculations and is 
divided in two steps. First, one has to get an approximate idea of the transition 
points using a direct extrapolation in 1/L of the numerical values of the gaps, 
computed at increasing L with a moderate number of DM RG states. Clearly, 
one may want to explore a rather large interval of values and so the increments 
in D will not be particularly small (0.1, say). Then, the analysis must be re- 
fined around the minima of the curves IS.E-YS-D with smaller increments in D 
and a larger value of M. In our problem, the approach that seems to give bet- 
ter results is standard finite-size scaling (FSS) theory [26, 29] (for instance as 
compared to the phcnomcnological rcnormalization group). 

Once the critical point is located, we take full advantage of the conformal 
structure by looking at the finite-size spectrum (see Eqs. (41) and (44)) of 
relevant and marginal operators. In practice, we select a number of states that 
tend to become degenerate with the ground state and look for straight lines 
in the AE-vs-L~^ plot. Then, from a best fit we expect to have a very small 
offset (ideally a zero gap in the thermodynamic limit) and a slope given by 
the scaling dimension d multiplied by the velocity prefactor, v, which is absent 
in the field-theoretical formulation but has to be determined (in terms of the 
microscopic parameters) in a lattice system. In the latter case, Eq. (41) should 
contain also a term CooL, Cqo being the energy density of the problem at hand. 
Actually, due to the prefactor v, we have to imagine a self-consistent procedure: 
Depending on the type of the transition wc have in mind (that is, depending 
on the central charge c), we stick on one or more levels in the spectrum that 
have exactly d= 1. Then the slope of these is nothing but v. Once the velocity 
is estimated, one uses Eq. (41) to best fit the product cv and see whether the 
value of c and the hypotesis on the universality class are self-consistent or not. 

To clarify the matter, let us start with the simpler case of the Haldane-Neel 
transition, that is thought to be in the 2D Ising universality class. Fixing A ~ 0.5 
we find £'c(0.5) = -1.2, and the /^-function method [29] yields i/(0.5) = 1.023 ± 
0.009, as far as the gap exponent, AE oc {D — DcY is concerned. Moreover, we 
observe the following nontrivial feature of the spectrum: The massless modes 
described by the CFT seem to be all and only the levels within S^^^ = 0, while 
those with S^^^ ^ mantain a finite energy gap in the limit of large L. Hence, 
the reference state for the calculation of v will be the second excited state in 
^tot — 0; corresponding to the primary field of conformal dimensions (1/2,1/2). 
Using quadratic extrapolations in 1/L we get v = 2.44, and consequently e^o = 
-2.0011961 ± 0.0000006 and c = 0.5008 ± 0.0008, thereby confirming the Ising 
universality class. The scaling dimensions can be estimated from the slopes of 
the straight lines in a plot like that of Fig. 4. In Table 1 the theoretical values 
anticipated in Sect. 4 are compared with these numerical estimates. The overall 
agreement is good (7 % in the worst case). Note that all the marginal operators 
have nonzero momentum and so they cannot represent a valid perturbation 
to the continuum Hamiltonian in as much as they would break traslational 
invariance. The absence of marginal operators suggests that each point of the 
Haldane-Neel transition corresponds to the same c = 1/2 theory and the line 
in the phase diagram is "generated" by the mapping from the discrete spin 
model to the continuum CFT. Repeating the same passages at A = 1 we get 
Dc{l) = -0.315, 1/(1) = 1.003 ± 0.006 together with v = 2.65, = -1.62651, 
c = 0.498 ± 0.002, that is, again a c = 1/2 contimmm theory. 

We now pass to an example of numerical investigation of a c = 1 line of 
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Figure 4: Energy differences, divided by 2it, plotted vs 1/L at the Ising tran- 
sition (A = 0.5,1? = —1.2). Points represent the numerical values obtained 
with multi-target DMRG runs (M = 216) collecting nine excited states within 
^tot — 0- Continuous lines are best-fit whose slopes are given in Table 1, to- 
gether with the theoretical predictions of the scaling dimensions (the labels on 
the right indicate the multiplicities, all correctly met). 



critical points, namely the transition from the Haldanc to the largc-_D phase. 
In the past [20], a similarity with the critical fan of the Ashkin- Teller model 
has been suggested. The operator content of this model arises from Ginsparg's 
orbifold construction [24] and consists of a number of _R'-indepcndent scaling 
dimensions plus the contributions coming from the pure Gaussian part (free 
boson) discussed in the previous Section. The fact that we do not observe K- 
indepcndent dimensions (apart from trivial secondaries of the identity) indicates 
that the continuum description of our spin-1 Hamiltonian with PBC at the 
Haldane-large-Z) transition should be purely Gaussian rather than "orbifold- 
like". 

In order to support this claim, we try again to match the whole spectrum of 
the relevant and marginal operators (d < 2). The difference with the c = 1/2 
case is that here wc have to fix not one but two nonuniversal parameters, v 
and K (see Eq. (50)). As regards the former, the velocity stems from the first 
and second excited states in S^^^ = 0. Note that in choosing these levels we 
are assuming, sclf-consistently, that i^T > 1 so that the two secondaries of the 
identity {d = 1) come first than the primaries with {m = 0,n = ±1), having 
do,±i = K. As far as the Luttinger parameter K is concerned, we have to inspect 
the spectrum in other sectors of Sf^^. too. In particular, the first excited state 
lies in \Sfgf \ = 1, that corresponds tom = ±l,n = Oin Eq. (50). The value of 
K is obtained from the slope d±i,o = 1/4:K in a plot similar to that of Fig. 4. 
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More generally, in order to check the self-consistency of the hypotesis c = 1, we 
have computed the finite-size spectrum of relevant and marginal operators in 
different sectors of S^^^ for a couple of critical points on the Haldane-large-£> 
line (first two rows of Table 2). Once that v and K arc numerically determined, 
the structure of the Gaussian spectrum is correctly reproduced (including the 
multiplicities) and the overall comparison is satisfactory since, in worst cases, 
the relative difference docs not exceed 3% (sec plots and tables of Rcf. [17]). 
The agreement with the theoretical predictions of the mapping in the planar 
regime is also remarkable. If we plug the coordinates of the critical points in 
the formulae of g and v for the Gaussian model derived above (Eq. (52)), we 
obtain v = g = 2.07, K = Tr/g= 1.52 at (A = 0.5, D = 065) &nd v = g = 2.45, 
K = Tr/g = 1.28 at (A = 1, £» = 0.99). 

Table 2: Velocity, central charge and ground state energy density for some 
critical points on the Haldane-large-D transition line. The numbers are the 
outcome of DMRG calculations with L = 16, 20, 24, 32, 48, 64 and M = 405, for 
cases with > 1, or L = 16, 20, 24, 28, 32, 36, 40 and M = 400, for cases with 
K < 1. The last two columns contain the numerical estimate of the nonuniversal 
parameter K, (according to the procedures described in the text) and the gap 
exponent obtained from the CFT formula v = 1/(2 — K). The error on Coo is 
typically of one unit in the last digit or better. 



[A,i?c(A)] 


V 


c 




K 


V 


(0.50,0.65) 


2.197 ± 0.004 


1.008 ± 0.003 


-0.908765 


1.580 ± 0.004 


2.38 


(1.00,0.99) 


2.588 ± 0.006 


0.997 ± 0.003 


-0.859152 


1.328 ± 0.004 


1.49 


(2.59,2.30) 


3.70 ± 0.04 


0.99 ± 0.01 


-0.675099 


0.85 ± 0.01 


0.870 


(3.20,2.90) 


4.445 ± 0.005 


1.133 ± 0.006 


-0.59132 


0.526 ± 0.007 


0.678 



Enforced by these quantitative predictions, we try to approach the multi- 
critical point where the c = 1 line meets the 0=1/2 one. Supposedly, the 
central charge at this point is c = 3/2 and it has been proposed [49] that the 
corresponding CFT is a SU(2)2 Wess-Zumino-Wittcn-Novikov model. If this 
was true, the two lines should join at the point where the effective Gaussian 
theory has K = \ [27] {FD point). Using A ^ Z? in the expression of g wc find 
that K = Ti I g(X) = 1 is satisfied for A ~ 2, while it is believed [13] that the 
multicritical point lies at A > 3. We guess that the two lines join at K <\, and 
in order to test this conjecture we study two more points: (A — 2.59, D — 2.30), 
again on the c = 1 line, and (A — 3.20, D — 2.90) proposed in [13] as the mul- 
ticritical point itself. Altough the steps are conceptually the same as above, 
here we encounter two additional complications. First, due to the closeness (or 
almost coincidence in the multicritical case) of the Ising transition, we observe 
the merging of the two (quasi)critical spectra. Hence, we have to target more 
states and separate the ones belonging to c = 1 from the ones belonging instead 
to c = 1/2. Second, we observe sizeable finite-size corrections from irrelevant 
operators. In fact, our analysis shows that we are moving at values of K smallcir 
than one towards K = 1/2 where certain irrelevant operators become marginal. 
As explained in [17], the last two rows of Table 2 are obtained by extracting 
K not from the first excited state, but rather from half the sum of the pair 
of levels with TO = 0,n = ±lin S^g^ = 0, to get rid of finite-size corrections. 
As anticipated, moving to the right of the Haldane-large-l? line the value of K 
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keeps on decreasing towards 1/2 {SD point) where we argue that this Hne meets 
the Haldane-Neel one and a first order transition starts. 

We close the section with a few comments on the hidden topological order 
measured by string order parameters (Eq. 23). It is expected that, leaving 
the Haldane phase, the Z2XZ2 symmetry is partially or totally restored. More 
precisely, when the c = 1 line is crossed, both Og and Og'^ vanish. As custom- 
ary, we can introduce two ofF-critical exponents that control the closure of these 
order parameters. For instance, fixing A and varying D about Dc{X): 

O's oc (Dc - Df^s , 0% oc (£)c - Df<^^ (54) 

Now, according to FSS arguments (sec. 5.1 of [24]), (3s and (3g arc related, 
via the gap exponent v, to their counterparts at criticality, that is, the scaling 
dimensions of the operators entering the associated string correlation functions. 
These dimensions, in turn, can be extracted from the slopes, rjs and in 
the log-log plots of 0'g^{D = Dc) evaluated at half of the chain. Using the 
relation 2ps = i^Vs i&nd analogously for the z channel) we find the values 
reported in Table 3 for a couple of critical points already discussed above. We 
should observe that the scaling dimensions r]s/'2 and ilg/2 are not contained 
in the c = 1 spectra cited above. However, we notice also that the numerical 
estimates of rjg are rather close to the values ^ '^"'"'"^ = K/2 and that these 
levels actually exist in the effective continuum theory provided that half-intger 
values of n are allowed in Eq. (50). In the XXZ spin-1/2 formulation this 
is known to correspond to twisted boimdary conditions on the chain. Thus, 
considering that the calculations presented here for the spin-1 case are with 
PBC, it's not surprising that the scaling dimensions associated with Og'^ are 
absent in the numerical spectra. Nonetheless, we believe that the closeness to 
K/2 is not accidental and in Ref. [17] we speculated about the possibility that 
the longitudinal string correlation functions (Eq. (21) with a = z) acquires, in 
the continuum limit, the asymptotic form 

g^{r) ~ (exp [TiVTTK^{0)] exp [±i\/^$(r)]) (55) 



so that the lattice string exp i'^J2]^=i^i 
uum twist operator exp [±WTTK^(r)]. 



is somehow related to the contin- 



Table 3: Exponents associated with the vanishing string order parameters at 
the Gaussian transitions taking place at the points indicated in the first column 
(see text for definitions). All the numbers are obtained with FSS on the data 
at L = 32, 48, 64. 80, 100 and M = 300. 



[A,i?c(A)] 


2(3s 


m 


Vs 


Vs 


(0.50,0.65) 


0.597 ± 0.009 


1.91 ± 0.02 


0.251 ± 0.002 


0.804 ± 0.003 


(1.00,0.99) 


0.407 ± 0.002 


1.10 ± 0.01 


0.2733 ± 0.0006 


0.741 ± 0.002 



6 Conclusions. 

In the present paper we have reviewed, to the best of our knowledge, part of 
the status-of the-art concerning Heisenberg spin chains, including biquadratic 
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interaction terms and various kinds of anisotropics, concentrating on the role 
of hidden symmetries in the various famiUes of spin models. We have discussed 
how the inclusion of anisotropy terms can drive the " standard" Heisenberg chain 
away from the Haldanc phase and how hidden symmetries (and their sponta- 
neous breaking) are of great help in classifying the "massive" (gapfuU) phases 
of the model. The location of the critical lines of the model has been accurately 
obtained numerically, confirming and extending earlier predictions [13]. 

The combined use proposed here of analytical (CFT) and numerical {DMRG) 
techniques to investigate the critical properties of the models has proved to be 
a rather successful strategy to clarify the nature and structure of the criti- 
cal phases of the models. Numerical simulation techniques (Monte Carlo and 
DMRG, to quote only the most known ones) are of more and more frequent 
and extended use in almost all branches of Theoretical Physics. A blind use of 
them can however be more dangerous than helpful in understanding the phys- 
ical properties of the systems for whose study they are employed. We believe 
instead that an "educated" use of numerical techniques in support of analyt- 
ical approaches, as described here, can result in a powerful synergy that can 
be of great help in understanding the physics of many problems in Theoretical 
Physics. 
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